From transcriptome to biological function: environmental stress in an ectothermic vertebrate, the coral reef fish Pomacentrus moluccensis

Background
Our understanding of the importance of transcriptional regulation for biological function is continuously improving. We still know, however, comparatively little about how environmentally induced stress affects gene expression in vertebrates, and the consistency of transcriptional stress responses to different types of environmental stress. In this study, we used a multi-stressor approach to identify components of a common stress response as well as components unique to different types of environmental stress. We exposed individuals of the coral reef fish Pomacentrus moluccensis to hypoxic, hyposmotic, cold and heat shock and measured the responses of approximately 16,000 genes in liver. We also compared winter and summer responses to heat shock to examine the capacity for such responses to vary with acclimation to different ambient temperatures.
Results
We identified a series of gene functions that were involved in all stress responses examined here, suggesting some common effects of stress on biological function. These common responses were achieved by the regulation of largely independent sets of genes; the responses of individual genes varied greatly across different stress types. In response to heat exposure over five days, a total of 324 gene loci were differentially expressed. Many heat-responsive genes had functions associated with protein turnover, metabolism, and the response to oxidative stress. We were also able to identify groups of co-regulated genes, the genes within which shared similar functions.
Conclusion
This is the first environmental genomic study to measure gene regulation in response to different environmental stressors in a natural population of a warm-adapted ectothermic vertebrate. We have shown that different types of environmental stress induce expression changes in genes with similar gene functions, but that the responses of individual genes vary between stress types. The functions of heat-responsive genes suggest that prolonged heat exposure leads to oxidative stress and protein damage, a challenge of the immune system, and the re-allocation of energy sources. This study hence offers insight into the effects of environmental stress on biological function and sheds light on the expected sensitivity of coral reef fishes to elevated temperatures in the future.

Transcriptome analysis of pigeon milk production - role of cornification and triglyceride synthesis genes

BACKGROUND : The pigeon crop is specially adapted to produce milk that is fed to newly hatched young. The process of pigeon milk production begins when the germinal cell layer of the crop rapidly proliferates in response to prolactin, which results in a mass of epithelial cells that are sloughed from the crop and regurgitated to the young. We proposed that the evolution of pigeon milk built upon the ability of avian keratinocytes to accumulate intracellular neutral lipids during the cornification of the epidermis. However, this cornification process in the pigeon crop has not been characterised. RESULTS: We identified the epidermal differentiation complex in the draft pigeon genome scaffold and found that, like the chicken, it contained beta-keratin genes. These beta-keratin genes can be classified, based on sequence similarity, into several clusters including feather, scale and claw keratins. The cornified cells of the pigeon crop express several cornification-associated genes including cornulin, S100-A9 and A16-like, transglutaminase 6-like and the pigeon 'lactating' crop-specific annexin cp35. Beta-keratins play an important role in 'lactating' crop, with several claw and scale keratins up-regulated. Additionally, transglutaminase 5 and differential splice variants of transglutaminase 4 are up-regulated along with S100-A10. CONCLUSIONS: This study of global gene expression in the crop has expanded our knowledge of pigeon milk production, in particular, the mechanism of cornification and lipid production. It is a highly specialised process that utilises the normal keratinocyte cellular processes to produce a targeted nutrient solution for the young at a very high turnover.

The invasive Northern Pacific seastar, Asterias amurensis, transcriptome

Here we provide the supplementary data collection for the upcoming publication of the larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. A. amurensis is ranked among the most potentially damaging invasive species in Australia and has recently expanded its range along the eastern mainland coast of Australia. As a first step to study the genetic basis of adaptive change and other important evolutionary processes during a contemporary invasive range expansion we de novo assembled and characterised the transcriptome.

Vibration mechanosignals superimposed to resistive exercise result in baseline skeletal muscle transcriptome profiles following chronic disuse in bed rest

Disuse-induced muscle atrophy is a major concern in aging, in neuromuscular diseases, post-traumatic injury and in microgravity life sciences affecting health and fitness also of crew members in spaceflight. By using a laboratory analogue to body unloading we perform for the first time global gene expression profiling joined to specific proteomic analysis to map molecular adaptations in disused (60 days of bed rest) human soleus muscle (CTR) and in response to a resistive exercise (RE) countermeasure protocol without and with superimposed vibration mechanosignals (RVE). Adopting Affymetrix GeneChip technology we identified 235 differently transcribed genes in the CTR group (end-vs. pre-bed rest). RE comprised 206 differentially expressed genes, whereas only 51 changed gene transcripts were found in RVE. Most gene transcription and proteomic changes were linked to various key metabolic pathways (glycolysis, oxidative phosphorylation, tricarboxylic acid (TCA) cycle, lipid metabolism) and to functional contractile structures. Gene expression profiling in bed rest identified a novel set of genes explicitly responsive to vibration mechanosignals in human soleus. This new finding highlights the efficacy of RVE protocol in reducing key signs of disuse maladaptation and atrophy, and to maintain a close-to-normal skeletal muscle quality outcome following chronic disuse in bed rest.

De novo assembly and characterization of the invasive Northern Pacific Seastar transcriptome

Invasive species are a major threat to global biodiversity but can also serve as valuable model systems to examine important evolutionary processes. While the ecological aspects of invasions have been well documented, the genetic basis of adaptive change during the invasion process has been hampered by a lack of genomic resources for the majority of invasive species. Here we report the first larval transcriptomic resource for the Northern Pacific Seastar, Asterias amurensis, an invasive marine predator in Australia. Approximately 117.5 million 100 base-pair (bp) paired-end reads were sequenced from a single RNA-Seq library from a pooled set of full-sibling A. amurensis bipinnaria larvae. We evaluated the efficacy of a pre-assembly error correction pipeline on subsequent de novo assembly. Error correction resulted in small but important improvements to the final assembly in terms of mapping statistics and core eukaryotic genes representation. The error-corrected de novo assembly resulted in 115,654 contigs after redundancy clustering. 41,667 assembled contigs were homologous to sequences from NCBI's non-redundant protein and UniProt databases. We assigned Gene Ontology, KEGG Orthology, Pfam protein domain terms and predicted protein-coding sequences to > 36,000 contigs. The final transcriptome dataset generated here provides functional information for 18,319 unique proteins, comprising at least 11,355 expressed genes. Furthermore, we identified 9,739 orthologs to P. miniata proteins, evaluated our annotation pipeline and generated a list of 150 candidate genes for responses to several environmental stressors that may be important for adaptation of A. amurensis in the invasive range. Our study has produced a large set of A. amurensis RNA contigs with functional annotations that can serve as a resource for future comparisons to other echinoderm transcriptomes and gene expression studies. Our data can be used to study the genetic basis of adaptive change and other important evolutionary processes during a successful invasion.

First comprehensive multi-tissue transcriptome of Cherax quadricarinatus (Decapoda: Parastacidae) reveals unexpected diversity of endogenous cellulase

The Australian freshwater crayfish species, Cherax quadricarinatus Von Martens, 1868, is an important commercial and invasive species that is also being increasingly used as a model organism to address important and interesting questions in crustacean biology. Through deep sequencing of the transcriptome of C. quadricarinatus from the hepatopancreas and four other tissues, we examine the evolution of endogenously transcribed cellulase genes and provide new insights into controversial issues regarding the nutritional biology of crayfishes. A cluster assembly approach yielded one of the highest quality transcriptome assemblies for a decapod crustacean to date. A total of 206,341,872 reads with an average read length of 80 bp were generated from sequencing the transcriptomes from the heart, kidney, hepatopancreas, nerve, and testis tissues. The assembled transcriptome contains a total of 44,525 transcripts. A total of 65 transcripts coding for carbohydrate-active enzymes (CAZy) were identified based on hidden Markov model (HMM), and a majority of them display high relative transcript abundance in the hepatopancreas tissue, supporting their role in nutrient digestion. Comprehensive phylogenetic analyses of proteins belonging to two main glycosyl hydrolase families (GH9 and GH5) suggest shared ancestry of C. quadricarinatus cellulases with other characterized crustacean cellulases. Our study significantly expands the number of known crustacean-derived CAZy-coding transcripts. More importantly, the surprising level of evolutionary diversification of these proteins in C. quadricarinatus suggests that these enzymes may have been of critical importance in the adaptation of freshwater crayfishes to new plant-based food sources as part of their successful invasion of freshwater systems from marine ancestors.

Monotreme lactation protein is highly expressed in monotreme milk and provides antimicrobial protection

Monotremes (platypus and echidna) are the descendants of the oldest ancestor of all extant mammals distinguished from other mammals by mode of reproduction. Monotremes lay eggs following a short gestation period and after an even briefer incubation period, altricial hatchlings are nourished over a long lactation period with milk secreted by nipple-less mammary patches located on the female's abdomen. Milk is the sole source of nutrition and immune protection for the developing young until weaning. Using transcriptome and mass spectrometry analysis of milk cells and milk proteins, respectively, a novel Monotreme Lactation Protein (MLP) was identified as a major secreted protein in milk. We show that platypus and short-beaked echidna MLP genes show significant homology and are unique to monotremes. The MLP transcript was shown to be expressed in a variety of tissues; however, highest expression was observed in milk cells and was expressed constitutively from early to late lactation. Analysis of recombinant MLP showed that it is an N-linked glycosylated protein and biophysical studies predicted that MLP is an amphipathic, α-helical protein, a typical feature of antimicrobial proteins. Functional analysis revealed MLP antibacterial activity against both opportunistic pathogenic Staphylococcus aureus and commensal Enterococcus faecalis bacteria but showed no effect on Escherichia coli, Pseudomonas aeruginosa, Staphylococcus epidermidis, and Salmonella enterica. Our data suggest that MLP is an evolutionarily ancient component of milk-mediated innate immunity absent in other mammals. We propose that MLP evolved specifically in the monotreme lineage supporting the evolution of lactation in these species to provide bacterial protection, at a time when mammals lacked nipples.

Epigenetic deregulation in pediatric acute lymphoblastic leukemia

Similar to most cancers, genome-wide DNA methylation profiles are commonly altered in pediatric acute lymphoblastic leukemia (ALL); however, recent observations highlight that a large portion of malignancy-associated DNA methylation alterations are not accompanied by related gene expression changes. By analyzing and integrating the methylome and transcriptome profiles of pediatric B-cell ALL cases and primary tissue controls, we report 325 genes hypermethylated and downregulated and 45 genes hypomethylated and upregulated in pediatric B-cell ALL, irrespective of subtype. Repressed cation channel subunits and cAMP signaling activators and transducers are overrepresented, potentially indicating a reduced cellular potential to receive and propagate apoptotic signals. Furthermore, we report specific DNA methylation alterations with concurrent gene expression changes within individual ALL subtypes. The ETV6-RUNX1 translocation was associated with downregulation of ASNS and upregulation of the EPO-receptor, while Hyperdiploid patients (> 50 chr) displayed upregulation of B-cell lymphoma (BCL) members and repression of PTPRG and FHIT. In combination, these data indicate genetically distinct B-cell ALL subtypes contain cooperative epimutations and genome-wide epigenetic deregulation is common across all B-cell ALL subtypes.

Urea-based osmoregulation in the developing embryo of oviparous cartilaginous fish (Callorhinchus milii) : contribution of the extraembryonic yolk sac during the early developmental period

Marine cartilaginous fish retain a high concentration of urea to maintain the plasma slightly hyperosmotic to the surrounding seawater. In adult fish, urea is produced by hepatic and extrahepatic ornithine urea cycles (OUCs). However, little is known about the urea retention mechanism in developing cartilaginous fish embryos. In order to address the question as to the mechanism of urea-based osmoregulation in developing embryos, the present study examined the gene expression profiles of OUC enzymes in oviparous holocephalan elephant fish (Callorhinchus milii) embryos. We found that the yolk sac membrane (YSM) makes an important contribution to the ureosmotic strategy of the early embryonic period. The expression of OUC enzyme genes was detectable in the embryonic body from at least stage 28, and increased markedly during development to hatching, which is most probably due to growth of the liver. During the early developmental period, however, the expression of OUC enzyme genes was not prominent in the embryonic body. Meanwhile, we found that the mRNA expression of OUC enzymes was detected in the extra-embryonic YSM; the mRNA expression of cmcpsIII in the YSM was much higher than that in the embryonic body during stages 28-31. Significant levels of enzyme activity and the existence of mitochondrial-type cmgs1 transcripts in the YSM supported the mRNA findings. We also found that the cmcpsIII transcript is localized in the vascularized inner layer of the YSM. Taken together, our findings demonstrate for the first time that the YSM is involved in urea-based osmoregulation during the early to mid phase of development in oviparous cartilaginous fish.

Transcriptomics-based analysis using RNA-Seq of the coconut (Cocos nucifera) leaf in response to yellow decline phytoplasma infection

Invasive phytoplasmas wreak havoc on coconut palms worldwide, leading to high loss of income, food insecurity and extreme poverty of farmers in producing countries. Phytoplasmas as strictly biotrophic insect-transmitted bacterial pathogens instigate distinct changes in developmental processes and defence responses of the infected plants and manipulate plants to their own advantage; however, little is known about the cellular and molecular mechanisms underlying host–phytoplasma interactions. Further, phytoplasma-mediated transcriptional alterations in coconut palm genes have not yet been identified. This study evaluated the whole transcriptome profiles of naturally infected leaves of Cocos nucifera ecotype Malayan Red Dwarf in response to yellow decline phytoplasma from group 16SrXIV, using RNA-Seq technique. Transcriptomics-based analysis reported here identified genes involved in coconut innate immunity. The number of down-regulated genes in response to phytoplasma infection exceeded the number of genes up-regulated. Of the 39,873 differentially expressed unigenes, 21,860 unigenes were suppressed and 18,013 were induced following infection. Comparative analysis revealed that genes associated with defence signalling against biotic stimuli were significantly overexpressed in phytoplasma-infected leaves versus healthy coconut leaves. Genes involving cell rescue and defence, cellular transport, oxidative stress, hormone stimulus and metabolism, photosynthesis reduction, transcription and biosynthesis of secondary metabolites were differentially represented. Our transcriptome analysis unveiled a core set of genes associated with defence of coconut in response to phytoplasma attack, although several novel defence response candidate genes with unknown function have also been identified. This study constitutes valuable sequence resource for uncovering the resistance genes and/or susceptibility genes which can be used as genetic tools in disease resistance breeding.

Analysis of human breast milk cells: gene expression profiles during pregnancy, lactation, involution, and mastitic infection

The molecular processes underlying human milk production and the effects of mastitic infection are largely unknown because of limitations in obtaining tissue samples. Determination of gene expression in normal lactating women would be a significant step toward understanding why some women display poor lactation outcomes. Here, we demonstrate the utility of RNA obtained directly from human milk cells to detect mammary epithelial cell (MEC)-specific gene expression. Milk cell RNA was collected from five time points (24 h prepartum during the colostrum period, midlactation, two involutions, and during a bout of mastitis) in addition to an involution series comprising three time points. Gene expression profiles were determined by use of human Affymetrix arrays. Milk cells collected during milk production showed that the most highly expressed genes were involved in milk synthesis (e.g., CEL, OLAH, FOLR1, BTN1A1, and ARG2), while milk cells collected during involution showed a significant downregulation of milk synthesis genes and activation of involution associated genes (e.g., STAT3, NF-kB, IRF5, and IRF7). Milk cells collected during mastitic infection revealed regulation of a unique set of genes specific to this disease state, while maintaining regulation of milk synthesis genes. Use of conventional epithelial cell markers was used to determine the population of MECs within each sample. This paper is the first to describe the milk cell transcriptome across the human lactation cycle and during mastitic infection, providing valuable insight into gene expression of the human mammary gland.

Spectral sensitivity of cone photoreceptors and opsin expression in two colour-divergent lineages of the lizard Ctenophorus decresii

Intraspecific differences in sensory perception are rarely reported but may occur when a species range extends across varying sensory environments, or there is coevolution between the sensory system and a varying signal. Examples in colour vision and colour signals are rare in terrestrial systems. The tawny dragon lizard Ctenophorus decresii is a promising candidate for such intraspecific variation, because the species comprises two geographically and genetically distinct lineages in which throat colour (a social signal used in intra- and inter-specific interactions) is locally adapted to the habitat and differs between lineages. Male lizards from the southern lineage have UV-blue throats, whereas males from the northern lineage are polymorphic with four discrete throat colours that all show minimal UV reflectance. Here, we determine the cone photoreceptor spectral sensitivities and opsin expression of the two lineages, to test whether they differ, particularly in the UV wavelengths. Using microspectrophotometry on retinal cone photoreceptors, we identified a long-wavelength-sensitive (LWS) visual pigment, a 'short' and 'long' medium-wavelength-sensitive (MWS) pigment and a short-wavelength-sensitive (SWS) pigment, all of which did not differ in λmax between lineages. Through transcriptome analysis of opsin genes we found that both lineages express four cone opsin genes, including the SWS1 opsin with peak sensitivity in the UV range, and that amino acid sequences did not differ between lineages with the exception of a single leucine to valine substitution in the RH2 opsin. Counts of yellow and transparent oil droplets associated with LWS+MWS and SWS+UVS cones, respectively, showed no difference in relative cone proportions between lineages. Therefore, contrary to predictions, we find no evidence of differences between lineages in single cone photoreceptor spectral sensitivity or opsin expression. However, we confirm the presence of four single cone classes, suggesting tetrachromacy in C. decresii, and we also provide the first evidence of UV sensitivity in agamid lizards.

Tripartite motif-containing 55 identified as functional candidate for spontaneous cardiac hypertrophy in the rat locus cardiac mass 22

BACKGROUND: Left ventricular (LV) hypertrophy is a risk factor for cardiovascular death, but the genetic factors determining LV size and predisposition to hypertrophy are not well understood. We have previously linked the quantitative trait locus cardiac mass 22 (Cm22) on chromosome 2 with cardiac hypertrophy independent of blood pressure in the spontaneously hypertensive rat. From an original cross of spontaneously hypertensive rat with F344 rats, we derived a normotensive polygenic model of spontaneous cardiac hypertrophy, the hypertrophic heart rat (HHR) and its control strain, the normal heart rat (NHR).

METHODS AND RESULTS: To identify the genes and molecular mechanisms underlying spontaneous LV hypertrophy we sequenced the HHR genome with special focus on quantitative trait locus Cm22. For correlative analyses of function, we measured global RNA transcripts in LV of neonatal HHR and NHR and 198 neonatal rats of an HHR × NHR F2 crossbred population. Only one gene within locus Cm22 was differentially expressed in the parental generation: tripartite motif-containing 55 (Trim55), with mRNA downregulation in HHR (P < 0.05) and reduced protein expression. Trim55 mRNA levels were negatively correlated with LV mass in the F2 cross (r = -0.16, P = 0.025). In exon nine of Trim55 in HHR, we found one missense mutation that functionally alters protein structure. This mutation was strongly associated with Trim55 mRNA expression in F2 rats (F = 10.35, P < 0.0001). Similarly, in humans, we found reduced Trim55 expression in hearts of subjects with idiopathic dilated cardiomyopathy.

CONCLUSION: Our study suggests that the Trim55 gene, located in Cm22, is a novel candidate gene for polygenic LV hypertrophy independent of blood pressure.